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Concept  Analysis  Corporation  (CAC)  markets  a 
general-purpose  CADCAM  system,  CADCEPT.  CAC  has  tailored 
custom  CADCEPT  features  for  Redstone,  using  the  CADCEPT 
Macro  language,  which  is  a  powerful  tool  for  geometric 
modelling  using  parametric  design  techniques.  This  system 
was  used  to  build  a  geometric  model  of  a  missile.  IVV' 
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Task  1  was  an  initial  problem  definition  task. 
Conversations  with  Mr.  Ed  VAughn  at  Redstone  resulted  in  the 
submission  of  a  sketch  of  a  typical  (unclassified)  missile 
and  specifications  for  input  format  to  Redstone's  SPIRITS 
analysis  program.  The  goal  was  defined  to  model  this 
missile  in  CADCEPT,  and  to  demonstrate  that  sections  could 
be  cut  and  fed  into  SPIRITS  in  the  proper  format.  It  was 
discovered  that  the  primary  reason  for  interest  in  CADCEPT 
was  its  ability  to  do  faceted  geometric  modelling,  and  that 
another  system  of  interest,  BRL  CAD,  from  the  Ballistics 
Research  Laboratories,  was  of  interest  to  Redstone 
personnel.  It  was  thought  at  the  time  that  BRL  CAD  might  be 
superior  to  CADCEPT  for  the  purpose  of  preparing  input  to 
SPIRITS.  CAC  has  since  acquired  and  implemented  a  copy  of 
BRL  CAD,  and  will  make  a  case  that  CADCEPT  will  be  a  nice 
complement  for  BRL  CAD,  especially  because  it  is  more 
interactive  than  BRL  CAD,  it  has  a  drafting  capability, 
interfaces  to  machining  postprocessors,  and  a  powerful  macro 
capability.  A  macro  is  somewhat  like  computer  software,  but 
is  really  the  product  of  an  end  user  with  an  application, 
written  in  the  CADCEPT  command  language.  The  medium  of 
representation  is  text,  and  the  instructions  are  interpreted 
(not  compiled)  by  the  system,  inline  with  other  user 
instructions.  The  main  difference  between  macro  writing  and 
other  steps  in  the  operation  of  CADCEPT  is  that  macros  can 
be  saved,  and  parameters  can  be  varied  easily  at  later 
dates.  In  addition,  macros  have  control  structures  similar 
to  structured  programming  languages.  Listings  of  macros  for 
modelling  the  sample  missile  and  for  other  purposes  are 
given  in  Figures  1,  7,  9,  11,  13,  15,  18,  22  and  26. 
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Task  2  was  the  preparation  of  demonstration  macros, 
which  illustrate  the  capability  to  handle  the  geometric 
shapes  and  design  problems  identified,  and  to  simulate  new 
software  required.  Information  necessary  to  commence  Task  2 
was  received  on  October  31,  1988.  It  was  decided  to  write 
macros  for  the  following  components: 
o  Hemispherical  nose 
o  Chlinder 

o  Transition  piece  from  cylinder  to  rectangular 
parallelipiped 
o  Rectangular  parallelipiped 

o  Group  of  airfoil  stabilizers,  modelled  as  symmetric 
Joukowski  profiles  in  cross-section, 
o  Holes  in  the  transition  piece. 
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It  has  turned  out  that  because  of  a  bug  in  the  UNIX 
version  of  CADCEPT  which  we  were  subsequently  forced  to 
use  because  of  the  failure  of  CAC's  VAX  computer,  and 
consequently  the  VAX  version,  we  have  not  been  able  to 
demonstrate  the  volume  subtraction  necessary  to  pierce  the 
holes  in  the  transition  piece. 

Task  two  has  been  further  elaborated  to  output  data  from 
CADCEPT  into  Redstone's  engineering  analysis  program 
SPIRITS,  and  information  on  this  requirement  was  also 
furnished  on  October  31.  A  program  has  been  written 
implemented  to  "thin"  section  data,  in  order  to  prepare  it 
for  input  to  SPIRITS.  A  conversation  with  Mr.  Vaughn  about 
the  subject  gave  information  that  sometimes  too  much  data  is 
produced  by  ray  tracing  and  other  programs.  CAC's  thinner 
will  preserve  the  character  of  sections,  while  descarding 
insignificant  data.  A  listing  for  the  FORTRAN  source  will 
be  found  in  Figure  24,  sample  input  file  in  figure  23,  and 
sample  output  file  in  figure  25.  Figure  22  shows  that 
CADCEPT  is  capable  of  interfacing  with  SPIRITS  in  the  sense 
that  it  can  produce  a  file  of  points  in  the  format  given  us 
by  Mr.  Vaughn's  October  31  submission.  Figure  26  shows 
further  the  reverse  process,  that  of  feeding  CADCEPT  with 
points  produced  by  other  programs.  We  might  add  that  this 
capability  is  extra  to  the  IGES  capability  which  CADCEPT 
also  has,  in  case  data  is  produced  not  by  just  any  arbitrary 
program,  but  by  a  CAD  system  which  also  has  IGES  capability. 

Results  of  Task  2  are  displayed  in  Figures  1-26.  It 
can  be  seen  from  perusal  of  the  listings  that  shapes  can  be 
changed  considerably  by  the  modification  of  a  few 
parameters,  thus  illustrating  parametric  design,  and 
modelling  without  drawings.  The  CADCEPT  database  is  more 
than  just  a  few  interrelated  files-many  components  are 
connected  in  logical  ways.  Macros  allow  cross  checking 
between  dimensions,  so  that  when  one  dimension  is  changed, 
others  can  be  made  to  change  automatically. 

Task  3  was  the  determination  of  software-hardware 
requirements.  The  extent  of  CAC's  involvement  in 
configuring  a  turnkey  system  and  being  in  a  consultative 
relationship  was  to  be  determined.  We  have  determined  that 
the  target  configuration  is  any  VAX  computer  running  ULTRIX, 
dec's  version  of  UNIX,  with  a  Tektronix  4207-compatible 
terminal.  Most  recent  information  from  Redstone  is  that  the 
4207 's  have  been  replaced  with  more  advanced  Tektronix 
terminals.  We  have  verified  that  these  are  upward 
compatible  with  the  Tektronix  4107,  a  slightly  less  powerful 
version  of  the  4207.  At  the  time  of  this  writing,  we  have 
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not  implemented  the  4107  support  code  for  UNIX.  CAC  will  do 
this  implementation  soon  on  CAC's  UNIX  machine,  for  CAC's 
own  benefit.  The  use  of  the  Tektronix  terminals  is 
straightforward  and  needs  no  special  support,  since  they 
will  operate  in  "4014  emulation"  mode.  The  major  benefit  of 
native  code  support  for  the  4107-4207  machines  as  far  as 
modeling  capability  of  CADCEPT  is  concerned  is  that  this 
enables  a  graphics  tablet  to  be  used.  We  understand  that 
Redstone  does  not  use  graphics  tablets.  Much  of  the  CAD 
world  seems  to  be  moving  in  the  direction  of  iconic  menus, 
which  require  only  the  terminal.  We  plan  to  implement  icons 
sometime  this  year. 

The  verification  of  whether  CADCEPT  will  run  on  the  VAX 
under  ULTRIX  is  being  performed  at  the  moment  of  this 
writing,  and  we  believe  that  any  problems  will  be  minor, 
since  as  we  understand  it,  DEC  ULTRIX  FORTRAN  is  more 
compatible  with  VAX-VMS  than  f77  FORTRAN  which  we  are 
currently  using.  This  is  because  VMS  and  ULTRIX  are 
manufactured  by  the  same  company.  Digital  Equipment 
Corporation.  The  UNIX  version  is  more  portable  than  the  VMS 
version,  because  we  have  eliminated  the  system's  dependency 
on  byte  ordering.  The  major  VAX-VMS-related  system 
dependencies  in  the  port  to  UNIX  were  changing  the  VMS  "QlO" 
functions  to  UNIX  system  calls,  which  required  interfacing 
FORTRAN  routines  to  UNIX  C-f unctions.  These  were  written  by 
the  Phase  I  Principal  Investigator  prior  to  submission  of 
the  initial  proposal  for  this  SBIR  contract. 

Also  as  part  of  Task  3,  we  have  verified  by  actual 
inspection  of  the  code  and  implementation,  as  well  as  by 
bringing  it  up  on  CAC's  in-house  UNIX  computer,  that  another 
CAD  package  of  interest  to  Redstone,  BRL  CAD,  is 
sufficiently  portable  and  well  engineered  to  run  on 
Redstone's  equipment.  We  have  included  the  knowledge  thus 
acquired  in  our  Phase  II  proposal. 

Task  4  was  the  preparation  of  system  recommendations. 
This  amounts  to  our  Phase  II  report.  In  summary,  our  system 
recommendations  are: 

(1)  BRL  CAD,  CADCEPT,  and  Redstone's  analysis  programs, 
such  as  SPIRITS,  should  be  interfaced  together  in  a 
tightly  coupled  fashion. 

(2)  Add  new  software  to  integrate  the  kinds  of  geometric 
entities  between  the  packages,  and  to  customize  them 
further  for  missile  geometry. 


3 


(3)  Provide  training  and  consulting  services  in  the  use  of 
the  new  system. 

(4)  Add  new  three-dimensional  display  devices. 

(5)  Implement  an  iconic  menusystem  throughout. 

(6)  Interface  to  the  Eagle  system,  a  large  integrated 
system  in  use  at  Eglin  Air  Force  Base  for  Computational 
Fluid  Dynamics.  CAC  has  acquired  a  copy  of  this  and  is 
developing  expertise  in  its  use.  CAC  personnel.  Dr. 
Glance  in  particular,  have  experience  in  Computational 
Fluid  Mechanics.  This,  combined  with  a  general 
knowledge  of  and  commercial  practice  in,  continuum 
mechanics,  show  that  CAC  is  well  qualified  to  implement 
the  use  of  Eagle  in  connection  with  all  the  rest  of  the 
software  mentioned  above. 

(7)  Interface  with  CAC's  structural  mechanics  programs,  in 
an  intelligent  and  specialized  manner. 

Task  5,  presentation  of  results  and  response  to  review, 
has  been  done  remotely,  by  the  installation  of  CADCEPT  on 
Redstone's  equipment  and  by  this  Final  Report  and  the  Phase 
II  proposal. 
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The  following  is  a  list  of  figures  documenting  the 
Phase  I  example  problem  results,  ''nascad''  means 
"“/nexus/nascad/nascad,  my“ .  In  the  following  instructions, 
commas  and  periods  after  CADCEPT  commands  are  for  these 
instructions  only  -  not  for  CADCEPT. 

Figure  1:  Listing  of  Macro  MISSILE.  This  is  the  main 
macro  which  invokes  other  macros  and  assembles  the 
missile.  Found  in  nascad/repairs/missile.com  .  cat 
it. 

Figure  2:  Isometric  view  of  the  result  of  running 
MISSILE.  Found  in  nascad/MISSILE.XYZ.  LOAD 
"MISSILE. XYZ'  and  SHOW  "ISO". 


Figure  3:  Zoomed  top  view  of  the  result  of  running 
MISSILE,  showing  the  Joukowski  profile  of  the  fins. 
While  data  for  Figure  2  is  LOADed,  do  SHOW  "T",  FIT, 
SHOW,  ZOOM  20  20  20  -10  -10  -1,  SHOW 


Figure  4;  Top  view  of  the  result  of  running  MISSILE. 
Figure  3,  FIT,  SHOW 


Figure  5:  Right  view  of  the  result  of  running  MISSILE. 
Figure  4,  SHOW  "R",  FIT,  SHOW 


Figure  6:  Front  view  of  the  result  of  running  MISSILE. 
Figure  5,  SHOW  "F",  FIT,  SHOW 


Figure  7:  Listing  of  macro  HEMI,  used  to  product  the 
nose  piece.  Found  in  nascad/hemi . com  .  cat  it. 


Figure  8:  Isometric  of  the  result  of  running  HEMI. 
Found  in  nascad/HEMI . XYZ  .  LOAD  HEMI.XYZ  .  LOAD 
HEMI.XYZ,  SHOW  "ISO",  FIT,  SHOW 


Figure  9;  Listing  of  macro  CYL,used  to  produce  the 
cylindrical  part.  Found  in  nascad/cyl . com  .  cat  it. 

Figure  10:  Isometric  of  the  result  of  running  CYL. 
Found  in  nascad/CYL.XYZ  .  LOAD  "CYL. XYZ",  SHOW  "ISO", 
FIT,  SHOW 


Figure  11:  Listing  of  macro  TRANPIECE,  used  to  produce 
the  transition  piece.  Found  in  nascad/tranpiece.com. 
cat  it. 


Figure  12:  Isometric  of  the  result  of  tunning 
TRANPIECE.  Found  in  nascad/TRANPIECE.XYZ.  LOAD 
"TRANPIECE. XYZ " ,  SHOW  "ISO",  FIT,  SHOW. 
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Figure  13:  Listing  of  macro  BOX,  used  to  produce  the  base 
component.  Found  in  nascad/box.com.  cat  it. 

Figure  14:  Isometric  of  the  result  of  running  BOX.  Found 
in  nascad/BOX.XYZ.  LOAD  'BOX.XYZ',  SHOW  'ISO',  FIT,  SHOW. 

Figure  15:  Listing  of  macro  FINS,  used  to  produce  the  set 
of  eight  fins.  Found  in  nascad/ repairs/ fins. com  and 
nascad/repairs/ jeffins.com.  cat  it. 

Figure  16:  Isometric  of  the  result  of  running  FINS  & 
JEFFINS.  Not  stored  as  complete  object.  Do  steps  in 
nascad/repairs/missile.com,  just  the  ones  for  the  fins  part, 
except  ATTACH  'MAIN'  instead  of  'MISSILE'. 

Figure  17:  Mathematical  analysis  for  Joukowski  profile. 

Figure  18:  Listing  of  macro  PROFILE,  for  Joukowski  profile. 
Found  in  nascad/profile. com.  cat  it 

Figure  19:  Typical  Joukowski  profile,  result  of  running 
PROFILE.  Found  in  nascad/PROFILE.XYZ.  LOAD  'PF4R0FILE.XYZ, 
FIT,  SHOW. 

Figure  20:  Sections  of  missile,  using  PIERCE  command.  Found 
in  nascad/repairs/MISSILE.PRC.  LOAD  'MISSILE. PRC' ,  SHOW 
'ISO',  FIT,  SHOW. 

Figure  21:  Missile  showing  section  planes.  Found  in 
/nascad/ repairs/MISSILE. CUT.  LOAD  'MISSILE. CUT' ,  SHOW 
'ISO',  FIT,  SHOW. 

Figure  22:  Listing  of  macro  OUTPNTS  showing  one  method  of 
getting  points  out  of  CADCEPT  into  a  file  in  SPIRITS  format. 
Found  in  nascad/ outpnts.com.  cat  it. 

Figure  23:  Listing  of  the  file  of  points  gotten  by  running 
macro  OUTPNTS  on  the  model  produced  by  PROFILE.  Found  in 
nascad/PROFILE.PNTS.  Cat  it. 

Figure  24:  Listing  of  thinner  program.  Found  in 
nascad/thin,  nascad/thin. f ,  nascad/thin.o.  A  standalone 
program,  which  prompts  for  input  and  output  and  thin 
tolerance. 

Figure  25:  Listing  of  set  of  points  produced  by  using 
thinner  program  on  file  shown  in  Figure  22.  Found  in 
nascad/ROFILE.PNTS.  cat  it. 

Figure  26:  Listing  of  macro  INPNTS,  showing  how  point 
files  can  be  read  into  CADCEPT.  Found  in  nascad/ inpnts.com. 
cat  it. 
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’  MAKES  A  MISSILE  FROM  COMPONENTS  STORED  IN'  COMPONENT  FILES 
'  GOTTEN  BY  LOADS  BELOW. 

'  FLAG  40  'ON'  I  HOLD  SCREEN  ON  TEXT  AREA  FELL. 

DVAR  ' ALL 


.  N  T  :  -iLH  M  1  OS  I  LE 
ADC  EL  I  '  0  0  0 

!''GG“'  ADDED  CELL  HEMT  ' 
■ico-.3 

ADTO  'CVL.DAr 
C/L 


AD CEL  I  'OYL  0  O  0 
M3G=  ■  ADDED  CELL  CVL ' 

MSG 

AD  1 0  '  TRANP I ECE . DAT ' 

TRANPIECE  5 
ATTACH  'MISSILE' 

ADCELI  'TRANPIECE'  O  0  0 
MSG='  ADDED  CELL  TRANPIECE' 
MSG 

ADIC  'BOX. OAT' 

BOX 

ATTACH  ' M I SS I LE ' 

ADCELI  'BOX'  0  0  0 
M3G=  '  ADDED  CELL  BOX ' 

MSG 

ADIC  'FINS.DAf' 

F I  NS 

ATTACH  'MISSILE' 

ADCELI  'FINS'  O’  0  0 
ADIC  'JEFFINS.DAT' 

JEFF I NS 

MSG^ ■  ADDED  CELL  FINS' 

MSG 

D'AAR  'ALL' 

STOR  'MISSILE. XYZ ' 

MSG='  ALL  DONE.  ' 

MSG 


0 

SDIC  MISSILE. DAT ' 

EXIT 

Ircjpaii'-sl  183) 

Figure  1 •  Listing  of  macro  MISSILE.  This  Is  the  main  macro  which  calls 
other  macros  to  assemble  the  missile. 


Figure  3«  Zoomed  top  view  of  the  rewult  of  running  MISSILE, 
showing  the  Joukowski  profile  of  the  fins. 


> 


Figure  4 


Top  view  of  the  result  of  running  MISSILE. 


Figure  5.  Right  Tiew  of  the  result  of  running  MISSILE 


running  MISSILE 
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TH  [  S  FILE  13  h effi i  .  c 

l-JILL  riAKE  HEMISPHERE  AT  CYLTO-' .  R 

INITIALIZE  PARAMETERS 

G  2L  ON'  !  MAKE  RESULTING 

I  NUMBER  SEGMENTS 
:  :.5  NUM3EG  !  NUMBER  SEGS  FOR 
2  .  00 1 


ADTU3  CVLRAD 


I  ZA" 


|M 


jp  PARAMETEF;S 


'  SHOULD  NOW  HAVE  THE  "DRIVE"  POLY  ON  LEVEL  0 

UNSA 
ADDS  1 

’  PTl  AND  PT2  WILL  BE  ENDS  OF  QUARTER  CIRCLE  ARC 
PT 1 1  -=CYLTOP  ( ,1  )  +CYLRAD 
PTir:=CYLTO?  (2) 

='T13-=CYLT0P  (3) 

PT21 ^CYLTQP ( 1 ) 

PT22=CYLTGP (2) 

P  '2 Y^CYLTOP  (3>  -‘-CYLRAD 
-  ■’■1==PT1  1  ,  PTl  2,  PTl  3 
••T2-rpT21  ,PT22,PT23 
ADDA  PTl  CYLTOP  PT2 

l-i/ 


0  10 


MAK 


C.  I 


HE  HEMI. 


aELA 
ADOS  0 

•"i  V  {■ ;cr  c 

VOLS 
..  MWaRD 
OVAR  'ALL' 

STGR  HEMI.XYZ' 

MSG='  NOW  EXIT,  RUN  INWARDPP.com,  GET  BACK  IN,  DO  LODVOL  ?■  SHOW 
!  MSG 


Q 

GDIC  HEM I. DAT' 
EXIT 

C.-ascadI  155) 


Figure  ?•  Listing  of  macro  HEMI,  used  to  produce  the  nose  piece 


Vi  C 

E  i::y_ETEPi 

TM[S  FILE  IE  cyl .com 
'  WILL  MAKE  RIGHT  CIRCULAR  CYLIl'iDER. 

'  INITIALIZE  PARAMETERS 

'  Z-  COORD  OF  BOTTOM  CF  OVLIMDER. 
MT^FO.  '  Z-  COORD  OF  TOP  OF  CYLINDER. 

C/LFAD-^-.  !  RADIUS  OF  CYLINDER.. 


ADDS  1 

ADDA  CyLBQT  CYLRAD 

CLEV  0  !  NOW  HAVE  LOWER  CIRCLE  ON  LEVEL  0 

COPY 

MOVE  TLONG 

CLEV  5  !  NOW  HAVE  UPPER  CIRCLE  ON  LEVEL  5 

S  E  L.  A 
ADDS  0 

RULE  0  5  10  '  MAKE  THE  CYLINDER. 

DELS  ■  GET  RID  OF  THE  CCMPCNENTS 

3ELA 

VOL  3 

INWARD 

DVAR  ' ALL ' 

STOR  'CYL.XVZ' 

'  MSG='  NOW  EXIT,  RUN  INWARDPP. com ,  GET  BACK  IN,  DO  LODVOL , 
MSG 


9DIC  ^ CYL. DAT' 
EXIT 

Inaccadl  15E) 


SHOW . 


Figure  9 


Listing  of  macro  CYL,  used  to  produce  the  cylindrical  piece 


Figiire  10,  Isometric  of  the  result  of  running  macro  CYL. 


i'1;''CF  D 

'  INITIALIZE  PARAMETEf^S 

HI  -EE  !  HEIGHT  QF  CIRCLE 

I  HEIGHT  QF  SQUARE 

I  RADIUS  QF  CIRCLE 

!  HALF-  DIAGONAL  QF  SQUARE 


CYLRrD-4. 


NUNDER  SEGMEN'I'S  IN  ALL  PQLY3CNS, 
NUNBER  GE3S  FOR  ;ARC  DISPLAYS. 

!  INCREMENT  FOR  LINE  SEGS  IN  POLY. 


CYl0OT=::‘,.O,H1 
=■'.  ^-cvlpad.o.h;, 

P2-‘.:  .  C -'LRAO,  HI 
P3-  -CyLRAD,!)„Hl 
P4-0  .  ■'  >  ,.RAD  ,  HI 
PT-FPDXRAD  .ChHC 
P:>-o .  FGXRAD  ,H2 
P -EG  X  PAD . O . H2 
FS-=0,  •RQXRAD ’h2 
i'iU'ISEG^TI 
SYS  ilQ  NUMSEG 
PALQNG=P5-P6 
PINC^PALCNG/NUMSEG  ! 
PEEG=PS 
SYS  2  .  ZjOOI 
DDC 

!  FORMAT  (aOAl,2X,I5) 


p 

£  TRANPIECE 

!  tq  make  MISSILE  TRANSITION  PIECE 
DUM--i:l 

TRN3TEF1  DUM 
DEFDRW  'TRANPIECE' 

"!-'OPUP  I  adds  straight  line  POLY  FROM  P5  TO  PS, 

!  WITH  NUMSEG  SEGMENTS 

SELA 

COPY 


ROTA 

COPY 

-90.  BOXTOP 

CYLBOT 

ROTA 

COPY 

-90.  BOXTOP 

CYLBOT 

ROTA 

SELA 

-90.  BOXTOP 

CYLBOT 

CGMNE 

;CT  P5 

DE!...S 

!  GET 

SELA 

!  WE  i 

CLEV 

JNSA 

0 

!  ON  1 

ADDS 

1 

SYS  1 

IS  4*F:UMSEG 

ADDA 

PL  CYL9QT 

'  E  ' 

CLEV 

SELA 

C* 

ADDS 

0 

RULE 

0  5  10 

DELS 

2,i~i  p. 

'v'Oi,..S 

INWARD 

M'.'C-  MOW  LXI'!',  RUN  I  NWAiRDFP ,  com ,  GET  BACK  IN,  ADIC  TRAMPIEri^ 
DVAR  ALL' 

STDR  -TPANP  I  ECE .  X  Y  2  ' 

'  MSG 


Figure  lla.  Listing  of  macro  TRANPIECE, 


Mccrj—  ' 


AMfN  nr-i 


f, 

•"RNSTEPl  DUM 
DEFDRW  ■PAR'^1' 

LQOVOL 

'jPFDRW  'HOLE' 

!  NOW  MAKE  A  TYPICAL  HOLE 

SYS  II.  a  4*NIJMSEG 

l-iCX^  •  CYLRAD+BQXRAD)  /2. 

HCENT=HCX ,0,H1 

HCENT!JP=MCX,0,H2 

HRAD= (BOXRAD-CYlRAD) /4. 

UNSA 
ADDS  1 

ADDA  HCENT  HRAD 

CLEV  1 1 

;...'NSA 

ADDA  HCENTUP  HRAD 
SLEV  12 
SF^LA 
ADDS  0 

RULE  !.  1  12  13 

DELS  :  GET  RID  OR  THE  COMPONENTS 

DVAR 

STDR  MISSILE. XYZ' 

SE  LA 

'v'QLS 

INWARD 

MSG=  MOW  EXIT  and  RUN  INWARDPP. com,  BACK  IN  AND  RUN  STEPS' 
'  TYPE  MSG  NUMSEG 

E 

E  STEPS 
CUN^  tl 
TRNSTEPl  DUM 

LOAD  'TRANPIECZ. XYZ •  FlgUT®  lib,  Macro  TRANPIBCE. 

SELA 


iia 


A'^TAC--'  =Ai=^l 
riDCELI  PCXTCP  '  -G'.E  ' 

rr  L.;i'"r  iir  ■ 

GTGF:  TRANP^ECE.  X /G  ' 

:'13G=  NOW  EXIT  AND  RUN  VQLUNEPP .  com  .  BACK 
'  type  NSG  NUMSEG 


IN  AND 


■'i  jiv)-;;:’  ! 

'■Ti:3  -EPl  DUN 

“j , r  T  Q  -  !' !  !z  C  r' V 


t 


3D  I  T 
EXIT 


AT 


Figure  lie.  Listing  of  macro  TRANPIECE,  concl. 


THIS  -TLE  IS  bc;i.com 

WILL  MAKE  RECTANGULAR  PARABLE- I FED. 
INITIALIZE  PARAMETERS 


H  1 

!  Z-- 

COORD  OF  BOTTOM  C~  BOX. 

1  7_ 

COORD  OF  TOP  OF  POX, 

TO  X R'l 

; 

HALF-  DIAGONAL  jF  BOX. 

L  L.  ■■= 

(■C)XPAD,0,H1 

‘  1  ’  L 

M=0  ,,0  ,H1 

r-,  1  • 

[•••'  '•■= ' }  V.  •  A. 

A  ..CN  3 

. ..  •.  ./  —  |— \  jc;i  t;  •  •— v  y  • 

‘"r]M 

j 

'T"!  +  A..nNG 

,  OOn<.  -\ 

L  E  r;  X 

ly  :  ■«  H 

O'-j'IA 

!'.■■■  T  ’f*'-  jin' 

W  PCX' 

ADDS 

L 

aDEL 

3CXRAD  0  HI 

0  BOX RAD  HI 

COPY 

ROTA 

-OQ.  BOTTOM 

BOX TOP 

COPY 

ROTA 

COPY 

. 90.  BOTTOM 

BOX TOP 

ROTA 

-90.  BOTTOM 

BOX TOP 

SELA 

ADDS 

0 

'  THE  ROTAs  DID  UNSA  ON 

CONNECT  STLO 

DELE  '  WE  NOW  HAVE  LOWER  SQUARE  POLYGON  ON  LEVEL  0 

SELA 

COPY 

MOVE  ALONG 


CLEV 

5 

!  WE  NOW  HAVE  UPPER  SQUARE  POLYGON 

SELA 

ADDS 

Q 

RULE 

0  5  10 

!  MAKE  THE  PRISM. 

DELS 

!  GET  RID  OF  THE  COMPONENTS 

SELA 

VGL3 

I NWAR 

D 

DVAR 

RiLL 

STOP 

' BOX. XYZ ' 

'  MSn=  NOW  EXIT,  RUN  INWARDPP.  com ,  GET  BACK  IN,  DO  ;..ODVOL,  SHOW. 
'  MSG 


SDIC  BOX. DAT 
EXIT 
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Figure  13*  Listing  of  macro  BOX,  used  to  produce  the  base  of  the  missile 
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Figure  14*  Isometric  of  the  result  of  running  BOX. 


E  -'TNES-^EPt 

r ;  -i  I  3  F  I.  L  h  i  o  ■;  i  ri  S  .  C  O  ;  T1 

^  wi:ll  iiake  a  missile  fin 

'  INITIALIZE  PARAMETERS 


PI 3FIN^0,0,30 

;  L  o  i  L."  L  'J  — ' ,  ^  n 
'w  J.  ,  -  ,  I"!  1 

EPfFN-0 . 0  ,  H2 
XOLT :!.  -10,0.0 


Z-PRIME  COORD  OF  OUTSIDE  OF  BIG  FINS. 
WIDTH 
THICKNESS 

'  CENTER  OF  SMALL  FINS 
CENTER  OF  BIG  FINS 
30  i  MINUS  EIGFIN 

I  -PRIME  POINT  jF  inside  OF  BIG  FINS. 

'  -'^■RIME  POINT  OF  OUTSIDE  OF  BIG  FINS 
'  AXES  FOR  ROTATION 


YGU'^L==W/2.  0,0,0 
VCLT2=W/2.0, 10,0 

■  L  ■  ^ 

SYS  IS  NUMSE3 


'  NUMBER  SEGMENTS  IN  ALL  POLYGONS. 

!  NUMBER  SEGS  FOR  ARC  DISPLAYS  PDLYGDNALI ZATION 


ZOEC-  SM ALF I N-B I GF I N 
DOWNAX-0, 0,-20 


E  FINS 

FINS3TEP1 

UNSA 

OEFDRW  '^INS' 

'  TO  MAKE  SYMMETRIC  JOLKOWSKI  PROFILE. 

S0=. 19245009  !  3aRT<3)/9.0  CADCEPT  DOESN'T  YET  HAVE  SORT. 

:=W/4.0  !  CIRCLE  RADIUS. 

EPS=Sa'»T/C  '  PROFILE  THICKNESS 

M=EPS*C 

A=M-!-C 

THETA=  0.0 

PD^ISO. 

rLPI==2. 0*PI 

TH I  NC-^TUP  I  /NUMSEG 

FNLPT  THETA 

P3E3=X ,Y,H1 

THETA=THETA+THINC 

'■NLPT  THETA 

PEND^^^X  ,  Y  ,H1 

PL  I  rL==.  001  ,  ,  001  ,  .001 

A DDL  PBEG  PEND 

F-'OXL^PEND-PLITL 

SOXR=:PEND+PLITL 

OILY  3CXL_  BDXR 

A'BES-F'EMD 

NTTMFS^MLMSEG-I 

FOR  •'  MW  1  1 ;  NT  I MES  >  FNOTHER 

SELA 

A  D  D  S  O 

CONNECT 

C;  TV  0  '  ON  LEVEL  O 

COPY 

MOVE  SPAN 


LL-I-V 

SELA 

ADDS 


WE  NOW  HAVE  LOWER  POLYGON  ON  LEVEL  5 

Figure  15a«  Listing  of  macro  FINS, 


iE._A 


£ 

E  FtJOTHER 

'  ADDS  ANOTHER  SEGNENT  "0  THE  POLYGON 
THETA^THETA+THINC 
FNUPT  THETA 


E  F^iUFT 

c  ■ 3  "F”  [:  □  9  (  r  r  i  f  t  a  ) 

33Ar3M"-  A*  A  CGG"^ 

C2-C  •■9 

— 5"'  '\j  f":  !'/* 

:■  N:^=A*i::ijGT . n 

X  =  PNU)*  (  ■  .  0  -i-C2D) 
F':v;w:=A*SIhj(THETA) 

''  -:FNW*  (  1 . 0-C2D) 

DVAR  TH  OENGM  PNN  COST  C2  C2D 


G 

SD  i;C  '  FINS,  DAT  ' 
EXIT 

CrsDcci.rsI  1S9) 


Fig  15b. 


Macro  FINS,  cont 


INS 


'C 


Mf 


^  i  L-  -v 


0  4 

'rF  INTO  CORRECT  WIND  POSITION 

V'l'  J  0  7.:0  0  0  0 

TE  "IN  TO  BE  PERP  TO  MISSILE  BOD/ 
?  !•  L  0  0  0  20  0  O 

l.iT  1  Q  Jtr  T‘ ‘Z*  ‘pi-jirr 


■  I  h  1 1"'  .•■••. 

l_> .  i- 1 
C  CT ! 

'  NOW 

•  thes 

MOVE 
COPY 
^  NOW 
MOVE 
SCAL 


r  :  'A  o 

! .  I  T  e  !  KH  I  1  r 

vv  1 1... L  f\'Jv: 


'HE  FINS  IN'^0  FORWARD  POSITION 


.  I  'J 

TO 


WILL  BE  THE  LARGE  ONES 


LARGE  FINS  WILL  BE  COPIED  ,  REDUCED  IN  SIZE 
('■'  0  —20 
.  7  O  0  1 0 


AMD  MOVED 


SDIC  'JEFFINS.DAT' 


L  r-epai  r  <s  I  190) 


Figure  15c 


Macro  FINS,  concl 


Figure  17.  Mathematical  analysis  for  Joukowski  profile 


As  ^goes  from  0  to  2/7",  point  (x,y)  given  as  follows  will 
trace  out  the  profile; 


x= 


a  cos<^  -  m  + 


cf_(a_co8<^_3_m) _ 

a^+"m^I"2am”coi^” 


y=  a  sim^- 


caj_sin^ _ 

a^”m'^-”2am’co8 


where. 


with  C  computed  from  tmax  in  figure,  s  - 


m=:£C,  and  a=  m  -»■  C  =  C(1+<£) 


b  '-■'WUFILb 

‘  TD  MAKE  SYMMETRIC  J0UKDW3KI  P^lOFILE. 

FORMAT  (FI  0.3) 

'^■ROMPT  ENTER  WIDTH,  F10.3'  W 
=^'RQMPT  'ENTER  THICKNESS,  F10.3’  T 

N!JMSFG  =  40  !  can't  make  this  very  i!;:..scn  bioger-  nc  roon. 

I  Note:  13  -format  dcssn't  wc-rl  (-f-or  prompt  -for  thi 

19245009  '  SGRT(3)/9.0 

'  CIRCLE  RADIUS. 


i_,  -W  '  4  .  (i 
EP:i=SDr 

o,;  -r'-i-r' 


ROFILE  THICKNESS 


■-H  I  NC-‘:‘UP  I  /  NUMSES 
^iijpT  A 

-  BEG=^X  ,  Y  .  ii' 
THETA^THETA+THINC 
NU-T  "HETA 
PEND^^X  ,Y,0 
PLITL=,.  001  ,  ,  001  ,  .  001 
ADDL.  PBEG  PEND 
BOXL=PEND~PLITL 
PQX9=:pC:NJ0  -fPL  I TL 
SELV  BO XL  BOXR 
P BE 3= PEND 
NTIMFS=NUMSeG-l 
FOP  ( MW 1 - 1 ; NT I MES )  NGTHER 
STOP  PROFILE. XYZ ' 

DVAR  'ALL' 

MSG=  '  ALL  DONE' 

MSG 


E  MOTHER 

ADDS  ANOTHER  SEGMENT  TO  THE  POLYGON 
THETA=THETA+THINC 
NUPT  THETA 


E  NUPT 
TH"T1 

CGST==CDS  (THETA) 
DENGM^-^::Ax-A-+M*M-2. 0*A-»-M*C0ST 

C2-C*C 

C2D=C2/DEN0M 
PMW.-:=a*COST-M 
X^^^PMW*  ( 1 . 0+C2D) 

P|v!W=A*S IN  (THETA) 

Y-:^PMW*  (  1 . 0-C2D) 

DVAR  TH  DENOM  PMW  COST  C2  C2D 


It 


3DIC  '  PROF T!._E.  DAT  ' 


EXIT 
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Figure  18. 


Listing  of  macro  PROFILEf  which  produces 
Joukowski  profiles. _ 


Figure  19.  Typical  Joukowski  profile,  produced  by  PROFILE 


Figure  20.  Sections  of  missile,  gotten  by  using  the  CADCEPT  command 
PIERCE,  In  general,  we  may  produce  input  for  SPIRITS  this  way. 


Figure  21,  Showing  locations  of  cutting  planes  for  sections  in  fig  20, 


'  T  f-i  I  :::i  "  1 1_  E!  IS  O  Ui  t  p  i  I  '1  S  .  C..  O  ffl 

^  WILL  PRODUCE  A  FORMATTED  FILE  -ROM  A  FROFILE. 
OPEL!  'PROFILE.  PMT3  ' 

FORMA-:"  vTFIO.T}  '  REMEMBER-  -ORMAT  WAMTS  ELANi- 
LlA  J  f- ROFILE.  XYZ  ' 


7 -WERE ‘ I > 

PijP  OUf’IPT 


STOP  '  r.L)-r=''  JT3.  XYZ  ' 

:OSG='  ALL  DONE  WITH  GUTEMl'S 
'  M  L  G 


E  DLM-^'- 
PT^rVpTX  I,  MW 
WRITE  PT 

E 

r. 

SDIC  'nurPMTS,.  DAT  ' 
E  :<  I T 
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Figure  22.  Listing  of  macro  OUTPNTS  showing  one  method  of  getting 
points  out  of  CADCEPT  and  into  SPIRITS. 


Ti  r.f'if'. 

2. 464 

0 .  C^C)‘0 

1  ■»  i  i  f  n  1 

2.  359 

0. 003 

C.  000 

2 , 1  87 

*;> .  1 1 

0 . 000 

1 . 755 

0.  024 

0. 000 

1 . 672 

0. 044 

0.  ooo 

1 . 346 

0. 070 

0 . 000 

0 . 937 

0 .  1  00 

0 . 0  0  0 

’  )  M  607 

0.  132 

0  -  0C>0 

.215 

0.  163 

L'f , 

0 , 179 

191 

0-  000 

0.564 

„  2  1  4 

0,  000 

0. 933 

'j  .  .2  I?  '6 

0,  000 

(•;i  ^  733 

0 «  000 

!  .  =^.’^0.1 

O  -  777 

.  coo 

1 . 965 

0 .  2  j.  0 

T) ,  ooo 

0 .  1 33 

0  -  OOC) 

2. 282 

0.  146 

0 . 000 

2.417 

0.  102 

C .  000 

2.498 

0.  052 

0 .  ooo 

2.526 

0  -  000 

0 .  ooo 

2. 498 

-0. 052 

0 .  ooo 

2.417 

-0. 102 

0 .  ooo 

2.282 

--0.  146 

0 . 0C)0 

2.097 

-0. 133 

0 .  ooo 

1 . 865 

-0. 210 

0 . 0C)0 

1 . 590 

-0.227 

0 .  ooo 

1 . 277 

-0 .233 

0.  ooo 

0 . 933 

-0.223 

0 .  ooo 

0.564 

-0.214 

0.  ooo 

0.  179 

-0. 191 

0 .  ooo 

•:) .  2 1 5 

-0. 163 

0 .  ooo 

0.607 

-0. 132 

(!) .  0  (!)<!) 

0.987 

-0. 100 

0.  ooo 

1 . 346 

-0. 070 

0.  ooo 

1 . 672 

-0. 044 

0 .  ooo 

1 , 955 

-0.024 

0 .  ooo 

2.187 

-0 ,  0 1  1 

0 .  ooo 

2,  359 

-0. 003 

0 .  (D'Od) 

2.464 

0.  OO'D 

'0 .  (!)<!)(!) 

2. 500 

0 .  OOd) 

'0 .  (!)<!)(!) 
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Figure  23.  Listing  of  the  file  of  points  gotten  by  running 
the  macro  OUIPMTS  on  the  file  produced  by  PROFILE. 


PROGRAM  USETHIN 
C 
C 

INTEGER  FLUN,  EFLUN 
REAL»9  X(3),  TOL,  XMAX 
REAL#8  POINTS (5, 128) 

CHARACTER*32  FNAME,EFNAME 

DATA  FNAME  /'  '/ 

DATA  EFNAME  /'  '/ 

EFLUN- 12 
FLUN- 11 
WR1TE(6,135) 

135  FORMAT! 'ENTER  INPUT  FILENAME:') 

READ (5, 136) FNAME 
WR1TE(6,137) 

137  FORMAT! 'ENTER  OUTPUT  FILENAME:') 

READ  !5, 136) EFNAME 

136  FORMAT !A30) 

OPEN ! UN I T-FLUN , F I LE-FNAME (1:30), FORM- ' FORMATTED ' , 

1  STATUS-' OLD ',ERR*30) 

GO  TO  32 

30  WRITE!6,138) 

138  FORMAT ('OPEN  INPUT  FILE  ERROR') 

GO  TO  999 

32  MPTS-0 

1 0  READ ! FLUN ,101, ERR-34 , END-998 )  X 

101  FORMAT (3F10. 3) 

MPTS-MPTS+1 
DO  77  MW-  1,  3 

POINTS (MW,MPTS)-X (MW) 

77  CONTINUE 

POINTS (4, MPTS)«0 
POINTS (5, MPTS)-0 
GO  TO  10 

998  CONTINUE 
TOL-  0.015 
XHAX-  100.0 

CALL  THIN(POINTS,MPTS,TOL,  XMAX) 

OPEN ( UN I T-EFLUN , F I LE-EFNAME ( 1 1 30  > , FORM- ' FORMATTED ' , 

1  STATUS- ' NEW ' , ERR-3 1 ) 

GO  TO  33 

31  WRITE (6, 139) 

139  FORMAT ('OPEN  OUTPUT  FILE  ERROR') 

GO  TO  999 

33  WRITE (EFLUN, 115)  MPTS 

115  FORMAT (14) 

DO  20  1-1, MPTS 

20  WRITE (EFLUN, 101)  (P0INTS(MW,1) ,  MW-1,  3) 

999  STOP 

34  WRITE (6,  667)  MPTS 

667  FORMATS  GOT  READ  ERROR.  MPTS-  ',  14) 

STOP 

END 

Figure  24a.  Listing  of  thinner  program.  A  standalone  program. 


SUBROUTINE  THIN(X,  MPTS,  TOL,  XMAX) 

C  THINS  5D  POINT  ARRAY  X  SO  THAT  CHORD  HEIGHT  IS  GREATER  THAN  TOL. 

C  USES  NO  TRIGONOMETRIC  FUNCTIONS  OR  SQUARE  ROOT. 

C  NPTS  IS  NUMBER  OF  POINTS  IN  ARRAY  X. 

C  THE  EXTRA  2D  IN  X  IS  AUGMENTATION  BY  U  AND  V  PARAMETERS. 

C  WILL  NOT  ALLOW  POINTS  TO  BE  DISCARDED  WHICH  ARE  CHORDALLY 

C  FARTHER  AWAY  FROM  LAST  POINT  THAN  XMAX.  THIS  ALMOST  MEANS 

C  THAT  THE  MAXIMUM  CHORDAL  POINT  SPACING  OF  THE  OUTPUT  IS  XMAX. 

C  X  AND  XMAX  ARE  REPLACED  ON  OUTPUT  BY  THE  THINNED  CURVE,  I.  E.: 

C  THINNING  IS  DONE  "IN  PLACE".  THIS  MAY  CAUSE  TROUBLE  ON  CONVERSION 

C  TO  ANOTHER  LANGUAGE  SUCH  AS  C,  WHERE  ARGUMENTS  ARE  PASSED  BY 

C  VALUE. 

C 

REALMS  X(5,128),  TOL,  XMAX 
INTEGER  MPTS 

REAL*8  Dl(3),  D2(3),  CH2,  CHORD,  TOLSQ,  Q1 ,  Q2,  D3,  DSQ 
INTEGER  NPTS,  L,  K,  LAST,  LI,  I,  M 

C  CH2  IS  CONSTANT  LENGTH  TOLERANCE,  SQUARED.  NOTE  THAT  WE  DO  ALL  LENGTH 

C  COMPARISONS  ON  SQUARED  LENGTHS,  SINCE  THE  SQUARE  FUNCTION  IS 

C  MONOTONIC.  WE  THUS  AVOID  COMPUTATION  OF  SQUARE  ROOTS. 

CH2-  XMAX  ♦  XMAX 
NPTS-  MPTS 

C  DONT  BOTHER  TO  THIN  LESS  THAN  11  POINTS. 

IF  (NPTS  .LE.  10)  GO  TO  999 
CHORD-  -1. 

TOLSQ-  TQL**2 

L-1 

K=2 

LAST-1 

C  PROGRAM  FROM  HERE  ON  USES  FOLLOWING  STRATEGY; 

C  INCREMENT  K  UNTIL  MAX  CHORD  HEIGHT  OF  RESULTING  POLYGON 

C  COMPOSED  OF  POINTS  L  THRU  K»1  EXCEEDS  TOLERANCE. 

C  THEN  WE  SAVE  THE  KTH  POINT  AND  START  NEW  POLYGON 

C  ALSO  AT  KTH  POINT. 

5  K-  K  +  1 

IF  (K.GT.NPTS)  GO  TO  200 
DSQ-  -1. 

Ql-  0. 

DO  6  M-1,  3 

DKM)-  X(M,  K+D-  X(M,  D 

6  Ql-  Ql  *  DKM)  «D1(M) 

C  DONT  SAVE  POINT  IF  VERY  CLOSE  TO  LAST  ONE. 

IF  (DABS(Ql)  .LT.  l.E-8)  GO  TO  5 
LI-  LAST  +1 

C  CHORD  IS  SQUARE  OF  CHORD  LENGTH  FROM  LAST  TO  K. 

CHORD-  (X(l,  K)  -  X(l,  LAST))#*2  + 

1  (X(2,  K)  -  X(2,  LAST))**  +  (X(3,  K)  -  X(3,  LAST))**2 

DO  10  I-  LI,  K 
Q2-0. 

03-  0. 

DO  7  M-1,  3 

D2(M)-  X(M,  I)  -  X<M,L) 


Figure  24b.  Thinner  program,  cont. 


7 


Q2-  Q2  +  D2(M)  *02 (M) 

Q3»  Q3  +  Dl(M)  *  D2(M) 

C  DSQ  IS  DISTANCE**  OF  X(*,  K)  FROM  CHORD  (X(*,  L) ,  X(*,  K  +  D), 

10  DSQ»  DMAXKDSQ,  02  -  03  *03  /  01) 

C  THROW  OUT  (DONT  SAVE)  POINT  K 

IF  ((DSO  .LT.  TOLSQ)  .AND.  (CHORD  .LT.  CH2) )  GO  TO  5 
C  SAVE  POINT  K 

L  *  L  +  1 
DO  30  N=  1,  5 

50  X(N,  L)=  X(N,  K  -  1) 

LAST=  K-1 
GO  TO  5 
200  L=  L  +  1 

DO  250  N»  1,  5 
250  X(N,  L)=  X(N,  NPTS) 

MPTS=  L 

999  RETURN 

END 


Figure  24c «  Thinner  program,  concl. 
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Figure  25.  Listing  of  points  produced  by  thinner  program  used  on 
file  of  points  shown  in  figure  22. 


E  'NFN'^S 

'  T:^l:s  E^ILE  IS  inpnts.cQfTi 

WILL  READ  A  FORMATTED  FILE  AS  A  POL 
WE  JUST  ASSUME  HERE  THAT  IT  HAS  AT 
'  THIS  DEMO  PROCEDURE  SHOULD  BE  FIXED 
OPEN  RCFILE.PNTS' 

FORMAT  (14)  '  REMEMBER-  FORMAT  WANT 


YGOM. 

LEAST  THREE  POINTS. 
UP  FGR  OTHERWISE. 


BL 


r“  I —  r 

;-»r  • 


READ  THE  FIRST  TWO  FTS,  MAKE  A 
NTLV  WE  HAVE  TO  CREATE  VARIABL 


0 


LINE,  SELECT  THE  SECOND 
S  FIRST  dEFCRE  WE  READ 


iD  NUj''P 


FORMAT  (3F10.3) 

":rAD  PTi 

A  p,  p.  j  IIT'  1“  ■*  p*  ~ 

PT  TNCO^^  .  001  ,  .  001  ,  .00: 

P  T !  E  T  =  =  P  T  2 — P T  N  C 

PTRIGHT=^PT2h-FTINC 

UNSA 

SELV  PTLEFT  PTRIGHT 

MAKE  UPPER  LIMIT  SOME  NUMBER  LARGER  THAN  LIKELY  NO  OF  P'^s 
'  IF  WE  MAKE  THIS  UPPER  LIMIT  VARIABLE,  CADCEPT  BOMBS. 

FOR  (MW=1 : NUMPTS)  MGREPT 


v'ERTE  V  . 
HEM . 


E 

E  MOREPT 
READ  PTI 

ADDV  PTI  '  ADD  THE  POINT  TO  THE  POLYGON  WHICH  WE  ARE  BUILDING 
IF  (RSTAT  -EQ.  'END')  MW=NUMPTS+1  !  GET  OUT  AT  END 
MSG='  IN  MGREPT.  MW,  PTt ,  RSTAT=  ' 

MSG 

MW 

PTI 

RSTAT 


sDIC  'INPNTS.DAT' 
EXIT 

Cnascad]  195) 


Figure  26 


Listing  of  macro  INPNTS,  showing  how  point  files  can  be  read 
into  CADCEPT. 


